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Abstract, This paper describes the methodology used by the General Rotorcraft 
Aeromechanical Stability Program (GRASP) to model the characteristics of the flow 
through a helicopter rotor in hovering or axial flight. Since the induced flow plays a sig- 
nificant role in determining the aeroelastic properties of rotorcraft, the computation of 
the induced flow is an important aspect of the program. Because of the combined finite- 
element/multibody methodology used as the basis for GRASP, the implementation of 
induced velocity calculations presented an unusual challenge to the developers. To 
preserve the modeling flexibility and generality of the code, it was necessary to depart 
from the traditional methods of computing the induced velocity. This is accomplished 
by calculating the actuator disc contributions to the rotor loads in a separate element 
called the air mass element, and then performing the calculations of the aerodynamic 
forces on individual blade elements within the aeroelastic beam element. 

Keywords. Aeroelasticity; Finite Elements; Helicopters; Induced Velocity; Rotary 
Wings. 


INTRODUCTION 

In September 1980, work began on developing the 
General Rotorcraft Aeromechanical Stability Program 
(GRASP). While numerous analyses (Ormiston and 
Hodges, 1972; Friedmann, 1973; Hodges, 1976, 1979; 
Warmbrodt and Friedmann, 1979; Friedmann and 
Straub, 1980; Davis et a/., 1974; Bielawa, 1976; John- 
son, 1977, 1980; Sivaneri and Chopra, 1982) are avail- 
able to perform aeroelastic analyses for rotorcraft, all of 
them are subject to major limitations (Johnson, 1986) 
in generality, flexibility, or theoretical consistency. The 
purpose for which GRASP has been developed is to 
provide a tool with enhanced capabilities that can be 
used to perform aeroelastic calculations for helicopters 
in hover and axial flight. 

The implementation of the hybrid finite-element/multi- 
body methodology (Hodges et a/., 1987a, 1987bb) in 


GRASP allows a structure to be modeled as a collec- 
tion of rigid bodies and flexible elements that can be 
connected in a completely arbitrary manner. While this 
methodology presents the analyst with a great deal of 
generality and flexibility in structural modeling, it also 
presents the developer with some challenges in imple- 
menting an appropriate representation of the helicopter 
flow field. Since the treatment of the flow around and 
through the rotor disk is an important part of any aero- 
elastic analysis of rotorcraft, it is vital that the induced 
velocity be calculated in a consistent manner. 

In this paper, current methods used to calculate the 
inflow will first be described. Then, the methodology 
used in GRASP will be discussed and the differences 
with the more traditional methods highlighted. Finally, 
the theoretical basis of the approach implemented in 
GRASP will be outlined. 


This paper is declared a work of the U. S. Government 
and is not subject to copyright protection in the United 
States. 
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METHODOLOGY 

Current aeroelastic stability analyses for helicopters use 
a variety of methods to calculate the steady-state and 


dynamic induced inflow. These range from simple, lin- 
ear models for uniform inflow to sophisticated, nonuni- 
form inflow models using free-wake analyses. While 
not breaking any new ground with respect to devel- 
oping new models, GRASP does take a different ap- 
proach with regard to its calculations of induced veloc- 
ity. Therefore, before describing the methodology used 
in GRASP, and its rationale, it will be instructive to 
look at some representative examples of the approaches 
taken in current analyses. 

Traditional Methods 

In many analyses (Ormiston and Hodges, 1972; Fried- 
mann, 1973; Hodges, 1976, 1979; Warmbrodt and Fried- 
mann, 1979; Friedmann and Straub, 1980), the steady- 
state induced velocity is calculated from a single, linear, 
closed-form expression that combines both momentum 
and blade-element contributions to the rotor forces and 
moments. This expression is a function of the rotor 
collective pitch angle (usually at the three-quarter ro- 
tor radius). Assuming uniform inflow, one takes the 
induced inflow velocity over the entire rotor disk to be 
constant with the same value as the theoretical value at 
the three-quarter rotor radius. Alternatively, one could 
assume that the inflow angle, which is the inflow veloc- 
ity divided by the local blade speed, is constant over 
the rotor radius with the same value as the theoretical 
value at the three-quarter rotor radius. 

Another approach to calculating the steady-state inflow 
is demonstrated by Davis et at. (1974), Bielawa (1976), 
Johnson (1977), and Sivaneri and Chopra (1982). In 
this method the inflow velocity is calculated as a func- 
tion of the thrust coefficient, which is usually given, 
however, the blade pitch angle required to produce the 
desired thrust is also a function of the inflow velocity. 
Thus, the computation of the induced velocity is nonlin- 
ear, and requires an iterative solution. The distribution 
of induced velocity over the rotor disk is then either as- 
sumed to be uniform, or specified by a set of assumed 
functions such as the Glauert induced velocity (Bielawa, 
1976). 

The method used in C AMR AD (Johnson, 1980) to cal- 
culate the induced velocity is more sophisticated than 
any of the preceding analyses. C AMR AD can use any 
of three methods to determine the induced flow. First, 
as above, a uniform inflow distribution is computed as 
nonlinear function of the thrust coefficient. Then, if de- 
sired, a nonlinear distribution can be determined from 
a prescribed- wake analysis, using the uniform inflow as 
an initial guess. If further refinement is needed, a free- 
wake analysis is performed using the prescribed- wake 
solution as the initial guess. 

Of the analyses just discussed, only a few (Johnson, 
1977, 1980) consider the effects of inflow dynamics. Ba- 
sic to this type of dynamic inflow analysis is the assump- 
tion that total forces on the rotor vary slowly enough 
that actuator disk theory is applicable to perturba- 
tion velocities. Comparisons with experimental results 
(Johnson, 1986) have shown that dynamic inflow can 
have a significant effect on aeroelastic phenomena. 

A common feature of all of the analyses discussed earlier 
is that the calculation of the steady-state inflow veloc- 


ities is performed separately from the main calculation 
of the steady-state deformation of the structure. Thus, 
the inflow generalized coordinates are not included in 
the state vector with the structural degrees of freedom. 
However, when inflow dynamics are included in the dy- 
namic problem (Johnson, 1977, 1980), those generalized 
coordinates are (and must be) included in the state vec- 
tor. Although this is a somewhat inconsistent treatment 
of the steady-state and dynamic inflow generalized co- 
ordinates, it does* not result in any significant analytical 
problems. This is a result of the coupling between the 
steady-state inflow and the structural deformations be- 
ing very weak. 

To separate the steady-state induced velocity calcula- 
tions from the structural calculations, contributions to 
the rotor loads from individual blade elements must be 
calculated at the same time as the flow field contribu- 
tions. Since the exact geometry of each blade may not 
be known, it is necessary to assume a relationship be- 
tween one or more blade parameters and the forces on 
the rotor. For example, in Sivaneri and Chopra (1982) 
the induced velocity is calculated as a function of thrust, 
and thrust is a function of the blade pitch angle at the 
three-quarter radius and the induced velocity. It is ap- 
parent in this case that there is some implied relation- 
ship between the blade pitch angle and the blade geom- 
etry and section aerodynamics. 

GRASP Methodology 

The axisymmetric flow field for a helicopter in hover 
or axial flight is represented in GRASP by an element 
called the air mass element. The inflow generalized co- 
ordinates associated with this element are then included 
in the steady-state model state vector as well as in the 
dynamic model state vector. This means the steady- 
state inflow velocity is calculated in parallel with the 
structural deformations, and that those velocities are 
fully coupled to the deformed state of the rotor blades. 
The inflow generalized coordinates are introduced into 
the model in a manner similar to that used to introduce 
structural degrees of freedom. That is, an air node is 
introduced to represent the flow field at a point on the 
axisymmetric axis of that flow field. 

One of the difficulties that arises from integrating the 
flow-field model into the structural model is related to 
the specification of the motion of the flow field rela- 
tive to the structure. It is known that if the rotor disk 
undergoes large deformations, the flow field will also un- 
dergo changes that follow and lag the disk deformations. 
This occurs because the flow field is not physically at- 
tached to the rotor, but is highly dependent (to say the 
least) on its location. However, since such large motions 
would result in periodic forces and moments for which 
GRASP will not account, that situation may be ignored. 
There are, therefore, two possible implementations that 
may be used. The first associates the flow field with 
an inertially fixed frame of reference. In this case, it is 
understood that large motions of the rotor (which are 
not allowed) will have no effect on the location of the 
flow field. The other option is to attach the flow field to 
the structure with the understanding that large defor- 
mations will violate the flow- field model assumptions. 


2 



For GRASP, the former option was chosen since it is 
closer to the actual physics of the phemonenon. 

Another difficulty with integrating the air mass element 
into the structural model arises because of the mul- 
tilevel substructuring capabilities, which enhance the 
flexibility and generality of GRASP in modeling com- 
plex structures. One of the concepts fundamental to the 
use of multilevel substructuring is that no substructure 
is required to have any specific knowledge of any sub- 
structures other than its parent. In the context of the 
flow-field calculations, the air mass element has no ac- 
cess to information on the geometry of the rotor. This 
makes it virtually impossible to make any assumptions 
that would allow the blade-element contributions to the 
inflow calculations to be included in the air mass ele- 
ment. Any assumptions that might be made would be to 
the detriment of the generality of the code. Therefore, 
the calculations of the momentum contributions from 
the actuator disk are separated from the blade-element 
calculations. The air mass element represents only the 
flow-field aerodynamics, while the blade-element aero- 
dynamics are isolated in the aeroelastic beam element. 

THEORETICAL DEVELOPMENT 

The theoretical development of the inflow equations is 
dependent on three components; the air node, the air 
mass element, and the aeroelastic beam element. The 
generalized coordinates that are used by GRASP to de- 
scribe the static state and dynamic perturbations of the 
induced velocities are supplied by the air node. The air 
mass element performs the calculations of the actuator- 
disk contributions to the inflow equations, while the 
aeroelastic beam element calculates the blade-element 
contributions. 

Under the assumptions used for this development, there 
are noncirculatory, blade-element contributions to the 
apparent-mass terms in the dynamic inflow. Some re- 
cent, but as yet unpublished work indicates that the 
dynamic inflow, apparent-mass terms result solely from 
circulatory effects. If this can be verified, some of the 
assumptions used in this analysis would have to be re- 
vised. 

Air Nod(? 

The induced velocity generalized coordinates are intro- 
duced into GRASP via the air node. These generalized 
coordinates are defined relative to an inertial frame of 
reference /, and they define the inertial air velocity at 
any point in the rotor flow field. Given that is an 
inertially fixed unit vector and A is also an inertial co- 
ordinate system with its origin at the center of the flow 
field (Fig. 1), the induced velocity U QI at a point Q 


U QI = -(ttf + ntr + (1) 

where r is the flow-field radial coordinate, and R is 

the position of Q relative to A in the 6, direction. 

7^, and 7(3 are the air node generalized coordi- 
nates. 


For the case of static inflow, generalized coordinates U* 
and 7j* r are used to represent uniform inflow velocity 
and radial velocity gradient at the center of the flaw 
field. The other two coordinates are not used. Dynamic 
inflow uses only generalized coordinates E/j 4 , 7^, and 
7(3 to represent the vertical and cyclic velocity pertur- 
bations. 

Air Mass Element 

The air mass element is implemented in GRASP to 
model the momentum flow of air through the disk of 
a helicopter rotor. In this element, the rotor is assumed 
to be an actuator disk, and the flow field a cylindrical 
region surrounding the disk (Fig. 1). The state vector 
for the air mass element is made up of the generalized 
coordinates for a single air node. In the following sub- 
sections the static and dynamic inflow models developed 
for the air mass element are discussed. 

Static Inflow. In the static case, the air is considered 
to be flowing steadily through the rotor disk. From 
momentum theory (Gessow and Myers, 1967), the dif- 
ferential thrust dT acting on a differential annulus of 
the rotor disk is related to the induced velocity U by 
the equation 

dT = 4 xp a rU |V + U\ dr (2) 

where p a is the air density, r is the rotor radial coordi- 
nate, and V is axial velocity of the rotor relative to still 
air (positive up). The total virtual work done by 
the thrust on the air is 

6W =4xp a J* U6P\V +U\rdr (3) 

where e is the root cutout radius, k is the rotor ra- 
dius, and SP is the virtual displacement of the air. The 
expression for virtual work is discretized by assuming 
that the induced velocity can be divided into a uniform 
velocity V* and a radial gradient 7* r so that 

U = Of+tf r r (4) 

The virtual displacement of the air is discretized iden- 
tically. Thus, 

SP = + s<t>* r r (5) 

When these expressions are substituted into the expres- 
sion for the virtual work, the coefficient of 6P? is equal 
to the rotor thrust while the coefficient of 6<f>^ r has the 
units of moment, but no real physical significance. 

Dynamic Inflow. The model for the inflow dynamics is 
taken from Pitt and Peters (1981). It is assumed that 
the freestream velocity of the rotor relative to still air 
is spatially and temporally uniform. This freestream 
velocity is augmented within the cylindrical region of 
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the flow field by the steady-state inflow velocity compo- 
nents just described. Then, infinitesimal dynamic per- 
turbations to the inflow are induced by dynamic per- 
turbations of the rotor thrust, roll moment, and pitch 
moment. 

The virtual work for the unsteady flow of air through 
the rotor disk is 


SMI = j* 2p a U\V + U\6Pdtl>dr 


HI ’- 6 


PaU6PdV ef t 


(6) 


K eff 


where V* is the rotor azimuth and is the effective 
volume of the cylindrical flow field. This statement of 
virtual work produces a system of first-order differential 
equations may be converted to a set of second-order 
equations and discretized by assuming 


U = Of +^iVr + f*i - $> l2 rsmil> + $ 13 rcosV> (7) 


SP = SP * — 6<f>* 2 r s\ntf> + 6</>* 3 rcosi/> (8) 

•A 

where Pj is the vertical perturbation of the induced 

velocity at the center of flow, ^> l2 and $ 13 are the cyclic 
perturbation gradients at the center of flow. 6P* is the 
vertical virtual displacement of the air at the center of 
flow, and and 6<j>* 3 are the cyclic virtual displace- 
ment components at the center of flow. 

Aeroelastic Beam E lement 

The aeroelastic beam element is the primary structural 
element in GRASP. It represents a slender beam that 
is subject to elastic, inertial, gravitational, and aerody- 
namic forces. Hodges (1985) derives the elastic, inertial, 
and gravitational forces in detail. This section will dis- 
cuss the derivation of the aerodynamic forces as they 
apply to the induced velocity calculations. 

In the following discussion, the symbol Q (for quarter 
chord) is used to denote the aerodynamic center. The 
static position of any quantity is identified as ( )', while 
( )" refers to the instantaneous position of the dynamic 
motions of the blade. As just mentioned, vectors are 
denoted by the underlined symbol. Measure numbers 
of vectors associated with a particular set of unit base 
vectors are subscripted with the identifier(s) for that set 
of unit base vectors. The unit base vectors used in the 
following discussion are shown in Fig. 2. 

The wind velocity vector W_‘ q at the aerodynamic cen- 
ter is calculated by subtracting the inertial structural 
velocity at Q" (V q r ) from the inertial air velocity at 
Q n ( If Q 1 ). In terms of the inflow generalized coor- 
dinates and V_ q 7 , the relative wind velocity measure 
numbers associated with the zero-lift-line basis vectors 
are 


wf»i = -cK‘ A (of + nt + K 

+ *&* &) - v$, 


( 9 ) 


The relative virtual displacement of an element of air 

qH 

with respect to the structure can then be ob- 

tained by applying Kirchhoff’s kinetic analogy to Eq. 9. 
All (’) quantities are replaced with 6( ), and all velocity, 
angular velocity, and velocity gradient symbols are re- 
placed by identically labelled virtual displacement, vir- 
tual rotation, and virtual displacment gradient symbols, 
respectively. All other terms are then discarded from 
Eq. 9. 

The magnitude of the relative wind velocity W at the 
aerodynamic center and the angle of attack a are time- 
dependent quantities that can be written in terms of the 
measure numbers of the relative wind velocity vector. 
Since this theory is two-dimensional, the relative wind 
velocity and angle of attack depend only on the measure 
numbers in the plane of the blade airfoil cross section; 
thus 


w = ^("#'1 y+(w&y (io) 

and 


tan a = 


W Q“ 
yy z >» 2 


(ii) 


o n 

The local airflow velocity gradient G z „ l2 I s also a time- 
dependent quantity that depends on the relative wind 
velocity. The subscripts 1 and 2 denote the gradient in 
the direction of the velocity measure number in the 
xi direction. This velocity gradient can be shown to be 


G q " - 
12 — 


dR% 


( 12 ) 


which, in terms of the inflow generalized coordinates 
and V_ q 7 , is 


n Q " - 

U Z” 12 ~ °n 


( 



dV Q,,/ 

aV Z"l 

(13) 


Like the virtual displacement, the virtual rotation of a 
structural element relative to the air 6T z „ 3 can be ob- 
tained by applying Kirchhoff’s kinetic analogy to Eq. 13. 
This is accomplished by replacing all (*) quantities with 
£( ), replacing all velocity, angular velocity, and flow 
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gradient symbols by identically labelled virtual displace- 
ment, virtual rotation, and virtual displacement gradi- 
ent symbols, respectively, and discarding all other terms 
from Eq. 13. 

Static Inflow. Since the relative wind velocity, the ve- 
locity gradient, and the angle of attack are all time- 
dependent quantities, contributions to the static part 
of the virtual work are obtained by separating out the 
static terms from Eqs. 9, 13, and 11. The static magni- 
tude of the relative wind is then 


where c is the local blade chord, W is the magnitude 
of the relative wind velocity, G^„ l2 is the flow veloc- 
ity gradient, and is the flow velocity normal to 

the zero-lift line (Fig. 2). The lift, drag, and moment 
coefficients (cj, c^, and c m ), respectively, are nonlinear 
functions of the blade angle of attack a. 

Dynamic Inflow. After the static quantities in Eqs. 9, 
11, and 13, have been removed, only the dynamic terms 
remain. The dynamic part of the magnitude of the rel- 
ative wind is written as 





(14) 


i a/Q* d/Q" wQ' i IrQ" 

Z_ w z'jW z h 2 

W 


(19) 


where the overbars indicate the static part. Similarly, 
the static value of the angle of attack is 


and the dynamic part of the angle of attack is 


tan d = 


rirQ' 

VLzx 

w Qt 


(15) 


& - 


u/Q* i foQ" _ i IrQ* v 1/Q” 

” Z'l™ Z n l ** Z'X™ Z"\ 

W 2 


(20) 


The expression for the virtual work 6W done by the 
aerodynamic forces over the length of the beam element 
is 


where the checks indicate the dynamic part. Then, the 
dynamic terms in the virtual work can be put into the 
form 


SV= f\-6S$“ i F z .. i + 6rf i ; a M)dx 3 (16) 
Jo 

where 6S% t{ is the virtual displacement of the structure 
relative to the air, 6T^„ S is the virtual rotation of the 
structure relative to the air, and Fz»% and M are the 
applied forces and moments at the aerodynamic center, 
respectively. The applied aerodynamic force vector F 
on a blade section (a distributed force per unit length 
of blade) is assumed to be 

£ = £ci\ + + Cnci I (17) 


where £ c is the circulatory lift, D is the drag, and £ nc 
is the noncirculatory lift. 

The equations that define the aerodynamic force compo- 
nents act on the aeroelastic beam element at Q and are 
determined from a quasi-steady adaptation of Green- 
berg’s thin-airfoil theory (Greenberg, 1947). 


Cc = \paW 2 CC, + 

V = ^p a W 3 cc d 

M = \paW 2 C*C m 



( 18 ) 



where 4$ &nd $s represent the displavement and ve- 
locity perturbations of all of the structural generalized 
coordinates. From this expression for the virtual work, 
the aerodynamic contributions to the aeroelastic beam 
element mass, damping, and stiffness matrices, M, C, 
and K , respectively, can be determined to be 
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M = AFH 

C = AEH + AFG + BH (22) 

K = AEG + BG + D 

Here M turns out to be symmetric, but neither C nor 
K are. Explicit expressions for the elements of M, C, 
and K can obviously be obtained by substitution. Such 
expressions are quite long and complicated: however, in 
view of GRASP’s method of evaluation of these matrices 
numerically from Gauss-Legendre quadrature, it is not 
necessary to obtain them. 

CONCLUDING REMARKS 

The method used in GRASP to model the rotor flow 
field for a helicopter has been described. The primary 
feature of this implementation that differentiates it with 
other approaches is the separation of the blade-element 
calculations from the actuator-disk calculations. Also, 
this method incorporates the inflow generalized coordi- 
nates in the state vector for the steady-state problem, 
which guarantees full coupling with the structural de- 
formations. 

Because of the approach used to implement the inflow 
calculations in GRASP, it is also possible to use im- 
proved flow-field models without having to develop an 
entirely new blade element. The analyst would then 
have at his disposal a prescribed- or free-wake flow field 
representation as in Johnson (1980), or perhaps an un- 
steady flow-field model like that recently developed by 
Peters and He (1987). However, since the current ver- 
sion of GRASP does not have an air node that is general 
enough to accomodate the different sets of generalized 
coordinates that would be required for these flow-field 
models, an improved, generalized air node would need 
to be developed. 
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